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ABSTRACT 

We present the first study of early dark energy cosmologies using N-body simulations 
to investigate the formation of non-linear structure. In contrast to expectations from 
semi-analytic approaches, we find that early dark energy does not imprint a unique 
signature on the statistics of non-linear structures. Investigating the non-linear power 
spectra and halo mass functions, we show that universal mass functions hold for early 
dark energy, making its presence difhcult to distinguish from ACDM. Since early dark 
energy biases the baryon acoustic oscillation scale, the lack of discriminating power is 
problematic. 

Key words: methods:N-body simulations — methods: numerical — dark matter — 
dark energy — large-scale structure of Universe 



1 INTRODUCTION 

Uncovering the nature of the dominant energy compo- 
nent of the Universe at the current epoch, dark energy, 
is a central goal of modern precision cosmolog y. The evi- 
dence fo r the existence of d ark en e rgy is strong (|Riess et al.l 
j 19981 ). iPerlmutter et al.l l|l999l ). iKomatsu et all l|2008l )). 
however on the theoretical front there are no lead- 
ing candidates despite a wealth of suggestions (see 
ICopeland. Sami fc Tsuiikawal (|2006l ) for a review). As ob- 
servational data becomes more precise and extensive due 
to current and future generations of large surveys and in- 
creasingly powerful instruments, there remains a number of 
theoretical challenges to determine with high accuracy the 
expected observational consequences of the myriad of possi- 
ble dark energy models. 

Dark energy affects the Universe primarily through 
the global expansion rate. This can be measured directly 
through dis tance measu r emen ts such as from Type l a 
supernovae jRiess et al. I l|l998t ). IPerlmutter et al.l (|l999t ). 
iKowalskT et al.l (l2008l vr over the last 10 billion years. How- 
ever, some dark energy models predict a non-negligible con- 
tribution to the early time expansion behavior, called early 
dark energy (EDE). To probe EDE requires observations 
tied to the very high redshift universe, such as the cos- 
mic microwave background (CMB) or baryon acoustic os- 
cillations (BAO). However, as shown in iLinder fc Robbersi 
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('2008'), early dark energy can be hidden in CMB observa- 
tions, even while shifting the intrinsic acoustic scale upon 
which BAO measurements rely upon. Thus failure to detect 
EDE can bias the cosmological model interpreted from CMB 
and BAO data. 

An alternative and complementary measure of dark en- 
ergy is through the observation of structure in the Universe. 
The change in the dynamical evolution of expansion relative 
to the ACDM model, say, alters the growth of structures 
in the inhomogeneous universe that we inhabit. The statis- 
tics of structure therefore encodes information about dark 
energy. Structure measures involve a different weighting of 
cosmological parameters than distance measures, thus the 
combination breaks degeneracies. Moreover, they can test 
the framework of gravitational growth of inhomogeneities, 
important for distinguishing between a physical dark en- 
ergy and changes to the laws of gravity. Finally, growth is 
a cumulative process, depending on the physical conditions 
from early times through the epoch at which the structure 
is observed. Thus, observations of structure are crucial to 
understanding cosmology at all epochs. Because of the im- 
plications of early dark energy for interpreting CMB and 
BAO measurements, and for understanding the fundamen- 
tal physics behind dark energy, the behavior of structure 
formation and evolution is a key point. 

Early dark energy models address the coincidence prob- 
lem by having the dark energy evolve relatively slowly com- 
pared to the dominant matter or radiation component; in 
some cases motivated by string theory they scale exactly 
with the background compo nent. Indeed, th is was one of 
the first dark energy mo dels |Wetterichlll988fl. ED E is dis- 
cussed in some detail in iDoran fc Robbersi l|20o3 ) and ref- 
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erences therein. Constraints from CMB and primordial nu- 
cleosynthesis data restrict the energy fraction in the early 
universe contributed by EDE to l ess than a few percent 
l|Doran. Robbers fc Wetterichll2007h but this could have sig- 
nificant effects. So me of th e impli cations for linear growth 
were examined in iLinded (|2006l ) and |d oran fc Robbers! 
ll2q06l) and extended semi-analytica l ly to non-linear theory 
bv lBartelmann. Doran fc WetterichI (|2006h . The main result 
of previous studie s of the non-linear growth in EDE cos - 
mol ogies (such as Bartelmann. Doran fc WetterichI (|2006| ) 
and iFedeh fc BartelmannI (|2007l )') is that for a fixed lin- 
ear theory matter perturbation amplitude at z = 0, (Tg, 
there are more collapsed structure in EDE models than 
ACDM . For higher redshifts this effec t increases, for in- 
stance iBartelmann. Doran fc WetterichI (|2006l ) found that 
at z = 1, there are up to an order of magnitude more dark 
matter halos for higher halo masses {Mhaio ^ 10^^ Mq). This 
can be simply interpreted as early dark energy, which does 
not appreciably cluster, suppressing the growth of structure 
at early times. Hence, to achieve the observed level of cur- 
rent structure, the early amplitude of clustering must have 
been greater. Other applications include the high Sunyaev- 
Zel'do vich amplitude poss i bly seen at high multipoles in the 
CMB ijSadeh et all l2007l : IWaizmann fc BartelmaiinI |2008|) 
and t he early formation of the first stars ( Sadeh fc Rephaelil 
l2008h . 

A key aim of this paper is to use N-body simulations 
to calculate the predicted number of dark matter halos in 
EDE cosmologies compared with ACDM. In particular, the 
non-linear level of structure is affected by the early growth 
in a way that can change predictions based on the usual lin- 
ear growth factor at a particular epoch. We examine both 
the halo mass function giving the abundance of collapsed 
objects and the non-linear matter power spectrum. N-body 
simulations are a robust method of determining these statis- 
tics and this paper critically examines predictions previously 
made analytically. 

In Section [2] we describe the N-body simulations we 
used to probe the effects of EDE on non-linear structure 
formation. We present the results for the abundance of dark 
matter halos in Section [3] and the results for the non-linear 
matter power spectrum in Sectionj?) We discuss the implica- 
tion of our results and make concluding remarks in Section 
111 



2 SIMULATION DETAILS 

The simulations were p erformed using the cosmological N- 
body code GADGET2 l|Springelll2005r ) suitably modified in 
order to model the expansion history for arbitrary dark en- 
ergy models. The initial CDM power sp ectrum was d eter- 
mined using the CMB code CMBEASY l|Doranll2005l ). Ini- 
tial realisations of the density field were created using part 
of the COSMICS code ijMa fc Bertschingerll 19951 ). modified 
to allow an arbitrary initial power spectrum. Given the real 
space displacement field, ^'o(x), generated from the linear 
power spectrum P{k) at redshift zero, the initial grid of 
particles are per turbed using the Zel'dovich approximation 
(|Zerdovichlll970l ) resulting in the positions and velocities at 
the starting scale factor of 
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where q are the real space grid co-ordinates and D{a) is the 
linear growth factor. The dots refer to derivatives with re- 
spect to cosmic time t. The linear growth factor and its 
derivative were calculated b y numerically inte g rating the 
growth equation detailed in iLinder fc Jenkins! l|2003l ). In 
all cases, the starting redshift was kept consistent between 
the different physical models. In general the linear power 
spectrum shape varies in the different models compared 
and therefore we made initial conditions separately for each 
model, always ensuring that the initial fluctuation amplitude 
returned the desired as at z = 0. The power spectra mea- 
sured from simulation outputs were compared to linear the- 
ory to ensure that the large scale power evolved as expected. 
All simulations and numerical analysis were performed on 
'The Green Machine' supercomputer located at the Centre 
for Astrophysics and Supercomputing, Swinburne Univer- 
sity. 

The early dark energy model is parametrised through 
its dimens ionlcss energy density as a function of redshift by 
(jDoran fc Robbers 2006) 
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in which and $1° are the density parameters of matter 
and dark energy today = 1 — in a flat universe 

as assumed here), wo is the equation of state of dark energy 
today and Sle is the energy density fraction of dark energy at 
early times, a <C 1. Note that our ile corresponds to n% from 
.Doran & Robbers (2006). Strictly for the initial comparison 
to (Bartelmann. Doran fc Wetteri^ (|2006l ). we also use the 
bending model of IWetterich deflned as 
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In both cases the dark energy equation of state is given by 
the stand ard formula w(a) = —(l/\3Qd(i — ^d )])dfld/dlna. 
Note that lDoran. Robbers fc WetterichI (120071 ) find that the 
marginal distributions of the likelihood for cosmological pa- 
rameters other than dark energy are insensitive to the spe- 
cific parametrisation. 

In all the simulations using the EDE form of Eq. |3]the 
cosmological model used was a flat universe with: Qrn = 0.3, 
fib = 0.048, h = 0.69, n = 0.98, erg = 0.76, and wo = -1, 
except where specifled otherwise. 

Dark matter halos in our simulation outputs were deter- 
mined using two different methods. The first is the Friends 
of Friends (EOF) method (jDavis et al.l!l985h usmg a con- 
stant linking length parameter of 6 = 0.2. An alternative ap- 
proach to defining halos in simulat ion outputs i s the spher- 
ical overdensity method (SO) (see iLukic et al.l l|2008l ) for a 
detailed comparison of EOF and SO halo methods). To find 
the SO defined halos, we employed the MPI version of the 
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AMIGA Halo Find e r B (A HF), successor of MHF intro- 
duced by iGill et all ||2C)04 ). In order to use this for EDE 
cosmologies, we made the appropriate modifications to the 
virial overdensity calculation. 

The matter power spectrum was calculated utilising the 
' chain ing the power' method, as described in ISmith et alJ 
using the cloud in a cell grid assignment scheme. 

The bulk of the simulations used 256'^ particles in a 256 
Mpc/h box starting at redshift z = 24. For each model we 
performed eight simulations using different realisations of 
the initial density field. The power spectrum and halo mass 
functions were averaged over all realisations. Additional sim- 
ulations were performed halving and doubling the box size 
as well as altering the starting redshift, softening length and 
accuracy parameters to ensure convergence. We also check 
that for fie approaching zero, the power spectrum and halo 
results converged to the results for constant w — wo ss ex- 
pected. The results of these tests indicated that the power 
spectrum ratio results are at most ~ 1% altered by different 
numerical parameters, with the most important factor being 
box size which affected the power spectrum ratios mainly at 
high k. Starting redshift and softening length make little 
difference to the ratios. The halo mass function results were 
somewhat affected by changing box s ize, o ur results agreed 
with the analysis of IPower fc Knebd (|2006l l for the effect of 
box size on the mass function. The ratio results fluctuated 
at the level of ~ 10% with varying box size. We do not claim 
state of the art precision in our mass function results, how- 
ever they are sufficiently accurate and converged to support 
our main thesis, as detailed in Section [S] 



3 HALO MASS FUNCTION 

The halo mass function (hereafter HMF), defined simply 
as the number of halos of a given mass expected per unit 
volume, n{M,z), is an important theoretical and observa- 
tional statistic of the dark matter density field. Indeed, the 
abundance is often claimed to be expone ntially sensitive to 
the co smology. The pioneering work of IPress fc Schecht"eil 
(| 19741 ) related the expected mass function to the variance, 
cr(_R) on some length scale R of the linearly evolved density 
field and the linear collapse parameter, Sc, which is derived 
from spherical collapse arguments and for ACDM is only a 
weak function of cosmology and redshift. 

The advent of larger N-body simulations allowed this 
mass function to be tested rigorously, resulting in improve- 
ments that modified the basic functional form to include 
ellipsoidal collap s e wit h parameters fit to simulation by 
ISheth fc TormenI (Il999l ). This in turn was i mproved upon 
by Jenkins et al. (|200ll ) and more recently bv I Warren et al.l 
l|2006h with an alternative approach of fitting directly the 
multiplicity function 
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by some function of a without reference to the collapse pa- 
rameter. Here p is the background matter density. In both 
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IJenkins et al] (120011 ) and lWarren et al.l l|200d) universal for- 
mulas for the multiplicity function, valid for a range of cos- 
mologies, were found. The effect of cosmology enters through 
the variance cr as a function of scale and redshift and the 
mean matter density. Writing the Press-Schechter style mass 
functions in this form gives results that are broadly in agree- 
ment with the Jenkins and Warren formulas with the latter 
being more accurate fits to numerical data, but qualitatively 
similar. Despite being formulated from ACDM, matter dom- 
inated and open cosmologies only, the Jenkins formula has 
been shown to be valid also for a large N-body simulation 
of a dynamical dark energy cosmology (|Linder fc Jenkini 
l2003h . 

In previous work on E DE cosmologies relying on estima- 
tion of the HMF, (such as [ Bartelmann. Doran fc WetterichI 
|2Q0dl) and lFedeli fc BartelmannI l|2007l )) the Sheth-Tormen 
mass function has been used to predict the rela- 
tive halo abundance of EDE and ACDM cosmologies. 
This required determination of the linear collapse pa- 
rameter for EDE cosmologies , for the derivation see 
iBartelmann. Doran fc WetterichI (|2006l ) . The analysis using 
this mass function suggested a significant increase in the 
number of massive halos in EDE models compared to ACDM 
models with the same as today. This was particularly true 
at 2 = 1. However, it is not clear that the Press-Schechter 
style mass functions, relying on the collapse parameter, is to 
be preferred for these cosmologies over the universal mass 
function approach of Jenkins and of Warren. For ACDM, 
the difference between the approaches is solely the degree 
of accuracy across a range of halo masses, but the functions 
broadly agree. 

However, in this work we compare the two approaches 
for EDE and find a marked qualitative difference in non- 
linear structure predictions, although they both use the 
same linear density behavior. We have reproduced Fig. 5 of 
iBartelmann. Doran fc WetterichI 12006) which showed the 
ratio of the HMF to the ACDM case for two EDE mod- 
els (b oth using the bending parametrisation of IWetterichI 
(|2004D . see our Eq. at z — and z — 1. However, in 
our reproduction (see Fig. [1} we have also included the 
prediction using instead the Jenkins and Warren formulas 
(note that these curves overlap and can barely be distin- 
guished on the scale shown). The Sheth-Tormen mass func- 
tions shown utilised the method for calculating the linear 
densi ty contrast at collapse time, Sc, f or ED E cosmologies 
from iBartelmann. Doran fc WetterichI l|2006l ). which results 
in a quite different value from ACDM. As can be seen, the 
predictions of the Sheth-Tormen and the Jenkins- Warren 
approaches are very different, and it is not clear a priori 
which approach should be used in the absence of rigorous 
calculation via N-body simulations. In this work, we fill the 
gap in the literature by performing these simulations and 
halo analyses. Note that the difference in the halo abun- 
dance ratios predicted by the Jenkins and Warren formulas 
is negligible for all cosmologies studied in this work and in 
subsequent figures we plot a single line that represents the 
prediction of both of these approaches. 



3.1 Press-Schechter vs. Jenkins- Warren 

As a direct comparison to the work of 

iBartelmann. Doran fc WetterichI l|2006h . we perform 
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iBartelmann. Doran fc Wettericl] l|2006h with the Ijenkins et al.l 

ood) 



1I2OOII ') and iwarren et alj (12009) mass functions also included. 



The Jenkins and Warren formulas are shown with dashed and 
dot-dashed lines respectively, however they are indistinguishable 
on this scale. They both clearly disagree with the Sheth-Tormen 
prediction, shown with dotted lines, for both models. For the 
Sheth-Tormen mass function, the upper curves are for z = 1 and 
the lower 2 = 0, however for the Jenkins- Warren mass functions 
the reverse is the case, the lower curves are for z = 1. 




M Mg/h 

Figure 2. Mass function ratios at z = 1 between EDE and 
ACDM for different methods of predicting the mass function. The 
solid lines show the simulation data with halos found with the 
EOF method. The dot-dashed lines show the SO method results 
for the same data. The dashed lines show the Jenkins mass func- 
tion (which gives almost the same prediction as the Warren for- 
mula). Finally, the dotted lines show the prediction of the Sheth- 
Tormen mass function. The Jenkins- Warren formula is clearly 
preferred by the simulation data. 



simulations of the same two EDE and one ACDM models 
they considered analytically. Model I has Q.m,oh^ = 0.146, 
VLbh^ = 0.026, h = 0.67, n = 1.05, mo = -0.93 and 
as = 0.82; Model II has iimfih"^ = 0.140, Qth^ = 0.023, 
h = 0.62, n = 0.99, wo = -0.99, and as = 0.78. Both 
models have an averaged value for the dark energy density 
during the matter dominated era of fld,sf = 0.04 (see 
iDoran et al.l l|200j l for the precise definition). 

The simulation results for the mass function, over- 
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Figure 3. Mass function ratios at z = for the models described 
in Section 13.11 The line styles are the same as Fig. [2] 



laid with the predictions from both the Sheth-Tormen and 
Jenkins- Warren mass functions, are shown in Figs. [2] and 
[3] The simulations show a good agreement with the predic- 
tions of the Jenkins formula and strong disagreement with 
the Sheth-Tormen function with the appropriate EDE linear 
collapse parameter. This has two crucial implications. The 
difference between ACDM and EDE in numbers of collapsed 
objects is much less than previously thought. Secondly, the 
abundance increase reverses at higher redshift, i.e. rather 
than an order of magnitude more halos in EDE at z = 1 as 
previously predicted, we find less halos at the high mass end 
for one of the models, and comparable number of halos at 
the low mass end. 

This result is important and troubling; instead of the 
number of halos being exponentially sensitive and a crucial 
tool for detecting EDE, it appears that it is quite insensitive 
to the existence of EDE at the few percent energy density 
level, at or below the upper limit satisfying current data. 
The implications are further discussed in Section [S] 



3.2 Early Dark Energy Mass Functions 

Although we used the bending model of EDE for dir ect com- 
parison to IBartelmann. Doran fc WetterichI (|2006l ) , an im- 
proved model of EDE w as formulated after that p aper. We 
treat the EDE model of iDoran fc RobbersI (|2006f ). given in 
Eq. |3l as the standard description of EDE in the remainder 
of the paper. It is somewhat more closely related to the par- 
ticle physics models and automatically goes to a constant 
energy density contribution fie at early times in both the 
matter and radiation eras. 

We compare the simulation results for this EDE model 
with the predictions of the Jenkins- Warren formulas in 
Figs. U and [S] Note that we do not plot the Sheth- 
Tormen predictions for the mass function, since we find 
that using the linear collapse parameter derivation from 



that using the linear collapse p arameter derive 
IBartelmann. Doran fc WetterichI ((2009 1 with the parametri- 
sation of Doran ( 20051 ) leads to a limit that does not converge 
and hence we cannot determine an appropriate 5c . In any 

^ If one assumes that the overd ensity of a spherical perturbation, 
A — 1, is proportional to a (as in lBartelmann. Doran &: WetterichI 
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Figure 4. Mass function ratios relative to ACDM at 2 = 1 for 
two EDE models of the form in Eq. |3]with different values of Sle 
but with the cosmology fixed otherwise. The dashed lines show 
the Jenkins- Warren mass function prediction for the ratios while 
the solid lines show the simulation results using the FOF method 
and the dot-dashed line the results using the SO method. 
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Figure 5. Mass function ratios at z = 0. The line styles are the 
same as Fig. |4] 



case, the results of the previous section indicated that this 
approach was not suitable for EDE cosmologies. 

As in Section fS.ll there are two key points to note, the 
first is that the difference between ACDM and EDE is neg- 
ligible at z = and at z = 1 (see Figs. |3] and [SJ, and second 
that the Jenkins- Warren mass functions correctly predict 
this result. Choosing to define the halos via the FOF or the 
SO method makes little difference to the ratio results. 

The simulations used in this study are insufficient to 
make a precision determination of the accuracy of the 
Jenkins- Warren formulas for EDE cosmologies, however it 
is clear that they are at least accurate to < 10% for predict- 
ing the relative mass function of EDE and ACDM and future 
work on EDE should use these rather than the linear theory 
collapse motivated mass functions. If required in the future. 



then since in the Doran-Robbers model D-|_(a) ~ 
^ then their ratio in Eq. 5 becomes a^^^c/^) and so 
5c formally goes to zero in the limit. To calculate the collapse 
parameter one actually needs the full numerical calculation. 



larger simulations with a greater mass resolution could be 
used to determine the mass function for EDE cosmologies at 
percent level accuracy, however this is outside of the scope 
of this work, and unnecessary at this stage due to the lack 
of precision in the corresponding observations. 



4 POWER SPECTRUM 

The mass power spectrum is central to extracting cosmo- 
logical information from galaxy redshift and weak lensing 
surveys. An important question is how does the addition of 
EDE change the non-linear mass power spectrum? Ideally 
we would like to be able, through simulations, to derive an 
understanding of the changes that EDE makes to the power 
spectrum, in order that we can predict this statistic for any 
EDE model without needing further lengthy simulations. 
The section presents the results of a systematic variation of 
cosmological parameters, in order to understand the non- 
linear power spectrum in EDE cosmologies. 

The approach we take is to examine how the ratio of the 
non-linear EDE power spectrum and ACDM compares to 
linear power spectrum ratios. In order to turn this informa- 
tion into a prediction of the EDE non-linear power spectrum 
in general requires accurate prior knowledge of the abso- 
lute value of the non-linear power spectrum for ACDM, for 
instance the simulation c a librat ed, halo model motivated, 
formula from lSmith et ahl (|2003t ) (know as Halofit), or any 
more accurate future method. In addition, EDE changes the 
shape of the linear CDM power spectrum. We therefore must 
calculate the correct lin ear power sp ectrum, using a pack- 
age such as CMBEASY (|Doranll200"5h . as used in this study. 
Armed with these two methods, we could then predict the 



non- linear power spectrum for EDE models if we could relate 
the ratio between the EDE and ACDM linear and non-linear 
power spectra. These will not generally be the same. This 
section addresses the ratio relation. 

Determining the absolute value of the power spectrum 
is a more intensive numerical t ask, outside the scope o f this 
work. However, as argued in [McDonald" et all (|2006l ) and 
iFrancis. Lewis fc Lindej l|2007h , many numerical errors will 
cancel in a ratio, and hence we can be confident, as long as 
convergence is demonstrated (see Section ^ that the ratio 
results are robust. Another issue is that the presence of EDE 
shifts the sound horizon and hence the l ocations in fc space 
of the BAO peaks in the power spectrum (|Doran et al.ll2007l : 
iLinder fc~R obbcrs 2008). However, our simulations concen- 
trate on the broad features of the power spectrum and do 
not have sufficient volume or statistics to resolve the percent 
level shift. Furthermore, the problem of how subtle non- 
linear effects alter the BAO pea ks in general i s an ongoing 
area of research (see for instance I Angulo et al.l |2008)). The 
focus of our current analysis is on the comparison between 
cosmologies of the broad features and small scales. 

A simple first step in predicting the non-linear power 
in EDE cosmologies might be to use the Halofit formula 
with the EDE linear power and growth factor. In the sub- 
sequent sections, we present the predictions of this modified 
Halofit form along with our simulation results to see how 
useful such an approach may be. In order to make this cal- 
culation the effective spectral index UcB and curvature C 
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Figure 6. Keeping erg constant while adding early dark energy re- 
sults in a remarkably similar non-linear power for k > 0.2 h/Mpc. 
This is despite the ratios in the linear power spectra, shown in the 
dashed lines, diverging. The modified Halofit prediction is shown 
by the dot-dashed lines. 
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Figure 7. At z = 1 models with early dark energy have slightly 
increased non-linear power on small scales, however the shape of 
the relative non linear power curve differs from the linear power. 

were determ i ned fr om the EDE linear power spectrum. See 
ISmith et al.l (|2003h for details and the Halofit formula. 

4.1 Varying Q.e 

We start with a standard concordance flat ACDM model 
with parameters as given in Section (2] We then consider 
two EDE models with different amounts of dark energy in 
the early universe, f2e, while keeping the rest of the cos- 
mology, including the growth normalisation today, ag , unal- 
tered. The linear and non-linear mass power spectrum ratios 
at z = and z = 1 are shown in Fig. |S]and[71 

As discussed in Section (2] the numerical convergence of 
the power spectrum ratios has been established at the ^ 1% 
level. The dependence of the ratio on numerical factors (pri- 
marily box size and particle resolution) is generally greater 
at higher k values. The results past k > 3/i/Mpc should 
perhaps be treated with more caution, but at lower k values 
the simulations are well converged. 

Compared to ACDM, for a fixed primordial spectral 



index n, the presence of dark energy in the early universe 
changes the higher k slope of the post recombination linear 
power spectrum, resulting in more large scale, but less small 
scale, linear power for the same overall normalisation ag. 
The growth amplitude in the early universe must be higher 
in EDE models than ACDM to compensate for the lower 
growth rate, in order to get the same erg today. The relative 
non-Hnear power between EDE and ACDM will therefore 
depend on those two governing effects, i.e. the change in the 
linear growth rate history and the change in the initial linear 
power spectrum shape. Non-linear power on small scales is 
coupled to the power on large scales and this non-linear 
amplification will be enhanced in the EDE models relative 
to ACDM due to the greater amplitude of large scale power 
(see Figs. |S] and [7}. The Halofit formula prediction (dot- 
dashed line in Figs. [6] and [Tj should be expected to take 
into account at least some of the non-linear amplification 
effects for the EDE linear power spectrum, since this formula 
was calibrated based on a variety of spectral shapes. While 
the correct linear growth factor at any given redshift can 
be calculated for EDE and used in the Halofit formula for 
the power spectrum at that redshift, the modified growth 
history of E DE models relative to ACDM, w ill n ot be taken 
into account. [Francis. Lewis fc Lindeil (|2003) and [Mi (|2007l ) 
found that the linear growth history in the early universe is 
an important factor in determining the non-linear growth at 
later epochs where the linear growth are matched. Therefore 
the scales where Halofit breaks down should grant an insight 
into the relative importance of the two factors, non-linear 
amplification and the linear growth history, as a function of 
scale. 

The results show that the prediction of the modified 
Halofit formula reproduces the simulation results in the 
translinear regime, around k = 0.3 — 1 /i/Mpc. As expected, 
the greater amplitude of large scale power in the EDE mod- 
els leads to enhanced non-linear amplification and hence the 
non-linear power ratio is greater than the linear. Beyond 
k > 1/i/Mpc, Halofit under-predicts the non-linear power, 
with the difference between this and simulations increas- 
ing with k. This can be understood from the ratio in linear 
growth history of the models shown in Fig. |S] The small- 
est scales are remembering the conditions of expansion, and 
linear growth, of earlier epochs, leading to a greater relative 
amount of non-linear growth at those scales. 

In order to decouple the effects of the modified growth 
history and spectral shape, we performed simulations sim- 
ulations in which the EDE models used the ACDM initial 
power spectrum while maintaining their correct expansion 
history. The z = result for these simulations are shown 
in Fig. |9] In this case Halofit predicts no difference between 
the models, due to the common linear power spectrum shape 
and normalisation. Examining Figs. |S]and|§]it can be seen 
that the difference between the simulations and Halofit pre- 
dictions in both figures are comparable. This reinforces the 
idea that Halofit correctly predicts much of the effects of the 
EDE spectral shape on intermediate scales and the modified 
growth history increases the power at small scales. 

This comparison is not perfect and we should not expect 
to be able to completely decompose the effects of spectral 
shape and linear growth history into a linear sum . These re- 
sults compare to those of lStab enau fc JainI (|2006D (see their 
Fig. 8) for the effects of a modified gravitational potential on 
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Figure 8. The linear growtli factors D{a) / D(l) of models with 
early dark energy relative to ACDM. The linear evolution of the 
growth is identical to within 1% after z ~ 1. 
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Figure 9. The same background cosmologies as Fig. E] but using 
the ACDM initial power spectrum for all models. The higher am- 
plitude of growth in the early universe in the EDE models leads 
to a greater non-linear power at small scales today. 



the non-linear power spectrum. In that case, the efltects on 
the non-linear growth of modifying the gravitational poten- 
tial were found to be negligible compared with the altered 
linear power spectrum shape. 

Note however that this is essentially the expectation 
previously proposed from analytic predictions of the halo 
mass function. The higher amplitude of growth in the early 
Universe in EDE models was expected to produce more col- 
lapsed structures at high redshift, with a residual increase 
in collapsed structures at low redshift due to the greater 
formation rate in the early Universe. Our power spectrum 
results do indicate an enhanced amount of non-linear growth 
relative to linear growth for the EDE models, however this 
enhancement is at the percent level rather than the con- 
siderable increase that would be seen if the relative halo 
abundance were an order of magnitude at high redshift as 
previously predicted. 

Interestingly, the differences in the non-linear growth 
rates due to the spectral shape and linear growth history 
appear to conspire to wash out the early dark energy sig- 
nal at z = 0. That is to say, the non-linear power in the 



EDE and ACDM models is matched to within a percent for 
scales smaller than k > 0.2/i/Mpc. Even ai z — 1 (Fig. [Tjl 
the high k power in the EDE models is well matched, to 
about (f2e/0.05)%. Since we expect fl^ < 0.03 from current 
observational constraints, this is a remarkable match. This 
result accords with the mass function results of Section [3.21 
that found essentially indistinguishable abundances of dark 
matter halos. Note that we are holding the cosmology apart 
from f2e fixed between EDE and ACDM. More generally we 
would not expect this match. For instance, a steep primor- 
dial spectral index n in an EDE model would compensate for 
the later flattening of the linear power spectrum due to the 
presence of EDE leading to a linear power spectrum after 
recombination much more similar to ACDM. In this case, 
as shown in Fig. we would expect a greater amount of 
non-linear power in the EDE model. 



4.2 Varying the Linear Growth History 

We now examine in more detail the influence of the lin- 
ear growth at different epochs on the non-linear power. In 
bottom-up hierarchical growth, the smallest scales reach the 
non-linear growth rate earlier than larger scales. We there- 
fore expect a relationship between the relative non-linear 
growth at different scales imprinted on the low redshift 
power spectrum and the alteration of the growth rate at high 
redshifts. This relationship was found to be important for 
scales with k > 1/i/Mpc in Section [4.11 To investigate this 
further, we have employed the ex tended EDE parametrisa - 
tion suggested in the appendix of lDoran fc Robber j l|2006l) 



+ ^2^(1 - a 



(6) 



where 7 controls the importance of the early dark energy 
terms at late times. Increasing 7 pushes to higher redshifts 
the epoch when this model changes from dark energy with a 
with nearly constant equation of state w = wo (i.e. ACDM 
when Wo = —1) to one with appreciable early dark energy 
density (or w approaching 0) . Th is extended parametrisat ion 
reduces to Eq.|3]when 7=1. In lDoran fc Robbers! l|2006l ') it 
was found that including the extra parameter does not alter 
the constraints on the other cosmological parameters when 
fitted to the data and hence this additional parameter was 
not deemed necessary. However, this parameter is useful in 
the current study as an aid to understanding the dependence 
of the relative non-linear growth as a function of the high 
redshift linear growth behaviour. The effect of varying the 
7 parameter on the linear growth factor is shown in Fig. [10] 

How does this altered linear growth behaviour change 
the non-linear growth? The models shown in Fig. [10] were 
simulated, with the measured non-linear power spectrum 
results shown in Figs. [11] and 1121 

At z = (Fig. I12p . models with a higher 7 parameter 
have less non-linear power at small scales than the model 
with 7 = 1. This can be understood from the relative linear 
growth histories of these models shown in Fig. 1101 Since the 
high 7 models remain more like the fiducial w = — 1 to higher 
redshift, the accumulated non-linear effects due to the linear 
growth history are smaller in this case than for 7 = 1. This 
causes the non-linear power to more closely track the Halofit 
prediction. 
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Figure 10. The effect of varying the 7 parameter in Eq. |6]on 
the hnear growth factor. Increasing the value shifts the departure 
from the fiducial ACDM to higher redshifts. 
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Figure 11. Power spectrum ratios with a varying 7 factor at 
2 = 1. The solid lines show the non-linear power while the dashed 
Unes show the linear power, both as a ratio to the ACDM model. 
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Figure 13. Power spectrum ratios relative to a fiducial too = 
—0.8 constant equation of state model. The solid lines show the 
non-linear power while the dashed lines show the linear power. 
Compared to Fig. [7] there is less difference between the linear 
and non-linear power spectrum ratios. 



At z = 1 (Fig. [11} Halofit matches the simulation results 
well for the high 7 models. Note that for both high 7 models, 
the linear growth factor evolution is essentially identical to 
ACDM in the period between z — 1 and z = (see Fig. 
nop . Despite this, the non-linear ratio of these models to 
ACDM does change between these two epochs. This increase 
is therefore due to the extra large scale power in the EDE 
models, and is predicted by Halofit. 

Altering the linear growth history through the 7 param- 
eter demonstrates the relationship between the high redshift 
linear growth history and the low redshift small scale power. 
As found in Section [4.11 this relationship, while important, 
is not strong enough to cause a significant increase in small 
scale power in EDE models. The small scale power at low 
and intermediate redshift is sensitive to the high redshift 
linear growth history at only the percent level. 
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Figure 12. Power spectrum ratios with a varying 7 factor at 
2 = 0. The solid lines show the non-linear power while the dashed 
lines show the linear power, both as a ratio to the ACDM model. 



4.3 Varying the Equation of State 

Until now we have compared EDE models to ACDM only. 
We will now investigate the effects from a different fiducial 
dark energy equation of state today, wo- We ran a series of 
models with wo ~ —0.8, with the rest of the cosmology the 
same as the ACDM model from Section 14.11 The relative 
power results showed very little dependence on wo, with 
the z = results almost indistinguishable from those with 
Wo — —1. At higher redshifts there is a weak dependence 
on wq in the non-linear power ratio, as seen by comparing 
the high k plateaus of Fig. [13] for wo — —0.8 vs. Fig. [7] for 
Wo = —1. 

In Fig. 1131 the ratios of the linear and non- linear power 
spectra are more similar than in Fig. [7] where wo ~ —1- This 
can be understood from the linear growth history which, 
compared to Fig. [S] shows a reduced difference in the linear 
growth history of the models at high redshift. This is simply 
due to the wo = —0.8 case dark energy acting more like 
matter in general (due to the less negative equation of state) 
hence the addition of the EDE component has less influence 
than in the ACDM case. 
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Figure 14. Power spectrum ratios relative to ACDM witli a 
higher (rg, 0.96, than that for Fig. |6] with erg = 0.76. With a 
greater amplitude of growth, non-linear effects are enhanced and 
the non-linear power spectrum ratios diverge more from the linear 
ratios. 



4.4 Varying the Growth AmpHtude 

The overall amplitude of structure is a major factor de- 
termining the rate of non-linear growth. It is therefore im- 
portant to understand how sensitive the relative non-linear 
growth is to this. We have performed a series of simulations 
varying the linear theory defined at z = 0. The result for 
ag = 0.96 (remembering the original model shown in Fig. 
[S] had as = 0.76) at z = is shown in Fig. [2] As might 
be expected, while ag, scales out in the linear power ratios 
(so the dashed lines are the same as in Fig. |6]), the higher 
amplitude of growth leads to enhanced non-linear effects. 
The difference in the non-linear and linear power spectrum 
ratios is increased with increasing erg. The results for z = 1 
have a similar increase in the power spectrum ratio and are 
not shown for brevity. 

Once again we can see that the modified Halofit for- 
mula matches the simulation results well at intermediate 
scales but under predicts the power in the EDE models on 
smaller scales. The scale where Halofit begins to under pre- 
dict the EDE power is pushed to a lower k value, again due 
to the increased non-linear effects from the higher growth 
amplitude normalisation. 

To test our intuition about effects governing the power 
spectrum results, we also ran simulations with fixed primor- 
dial power spectrum amplitude As rather than ag. The z = Q 
result is shown in Fig. 1151 For these simulations, one gets dif- 
ferent linear power today (and so different erg) but the early 
epochs of collapse, showing up in the very nonlinear (higher 
fc) regime today, are more similar. This can be seen in the 
power for A; > 1 /i/Mpc in the EDE models which increases 
with fc, heading back towards the ACDM power since the 
amplitude of power in the different models was more similar 
when those scales collapsed than it is today. 

Matching erg appears to be sufficient to ensure that the 
small scale non-linear power between ACDM and EDE mod- 
els are matched to within a few percent, regardless of the 
amount of early dark energy, Q.^ (see Fig. [6]) or the redshift 
when this EDE term turns on (see Fig. I12|l . This match is 
weakly dependent on erg; however varying erg by 0.2, a con- 
siderable amount, shifts the ratio in power by only a percent. 
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Figure 15. Simulations with fixed primordial amplitude As 
rather than fixed erg today. 

These results leads to a useful tool for confronting EDE 
models with power spectrum measurements, such as galaxy 
redshift surveys or weak lensing. The Halofit formula, modi- 
fied by using the EDE linear power spectrum and growth fac- 
tor, predicts the non-linear power ratio to within a percent 
out to fc ~ 1 ft/Mpc. For smaller scales, using the ACDM 
non-linear power matches the EDE result to within a few 
percent or so, particularly at low and zero redshift. This 
prescription is fairly loose but could be improved upon with 
a more quantitative study of the variation in power with 
cosmological parameters. The main point from this study 
is that the variations between EDE and ACDM are at the 
few percent level only, rather than there being a significant 
difference in non-linear power. 



4.5 Growth vs. Geometry 

The expansion history of an EDE model can be fit well out 
to 2 = 2 with a time varying equation of state w{a) = 
Wo -\- Wail — a) with Wa ~ 5r2e. Fig. 1161 shows the dark en- 
ergy density and equation of state for an EDE model and 
the matching {wo,Wa) model. Distance measurements, such 
as from Type la supernovae, will agree in the two models 
to 0.02% out to z = 2. However, the early universe histo- 
ries were very different. Does the growth history distinguish 
between these two models? To investigate this question, we 
simulated a pair of models, one using the Eq.[3]form of EDE 
with fie = 0.03 and wo — —0.95 and one using the {wo,Wa) 
form of dark energy with wo — —0.95 and Wa = 0.15 with 
the cosmology matched otherwise. The power spectrum ra- 
tio results are shown in Figs. 1171 and [181 

In this case the EDE model is taken as the fiducial, 
and we display the ratio of the {wo,Wa) model to this. Once 
again, we see that the non-linear power on small scales is well 
matched between the two models, despite the linear power 
ratios not matching at these scales. At z = 1 the difference 
between the linear and non-linear ratios is less, with more 
non- linear power at small scales in the {wo,Wa) model. Once 
again it is instructive to examine the linear growth history, 
shown in Fig. [19] 

This figure shows that the EDE model has a higher 
amplitude of growth in the early universe than the {wo, Wa) 
model. This is due to the greater amount of dark energy 




Figure 16. Dark energy density and equation of state comparison 
for an EDE and (wo,Wa) model with distance history matched 
to 0.02% out to z = 2. 
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Figure 17. Ratio results at z = for the comparison between 
the two parameterisations for dark energy given by to (a) = uiq + 
WaiX — a) and Eq.|3l where the parameter values are chosen match 
to the distance history. 
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Figure 18. As Fig.[T3but at z = 1. 



present at this time, leading to less linear growth over the 
whole history of the universe, requiring the EDE model to 
start at a higher growth amplitude in order to obtain the 
same linear theory crs today. The linear growth factor is very 
well matched between the two models after z ~2 (a~ 0.33), 
yet as can be seen, the non-linear growth between z = \ 
and z = is different in the two models, again due to the 
nonlinearity amplification effect. 



5 CONCLUSIONS 

In this work we have examined the effects of early dark en- 
ergy on the non-linear growth of structure through N-body 
simulations. The key result is that non-linear large scale 
structure, including the abundance of clusters, is relatively 
insensitive to the presence of EDE. This result is in contrast 
to previous analytic work which expected a substantial in- 
crease in non-linear structure compared to ACDM with the 
same erg. 

The implications of the results presented here are 
that EDE is much more difficult to detect than previ- 
ously thought. This presents a problem for dark energy 
cosmology, not only for the sak e of understanding dark 
ener gy but because, as shown in iDoran et al.l (|2007h and 
.Linder fc Robbers! (|2008l ). the presence of EDE can signif- 
icantly bias distance measurements using baryon acoustic 
oscillations if the possibility of EDE is not considered. Fit- 
ting for EDE however significantly reduces the discriminat- 
ing power of BAOs. 

Ideally it was hoped to det ect the presence of 
EDE through cluster abundances (see lFedeli fc BartelmannI 
(20071)), such as with future Sunyaev-Zel'dovich, weak lens- 
ing. X-ray, or optical surveys. This would provide an inde- 
pendent measure of EDE thus restoring the discriminating 
power of BAO measurements. However, this work demon- 
strates that this is not the case. 

As the amount of EDE increases, it becomes easier 
to distinguish in that the primordial amplitude of mat- 
ter density perturbations must be increased to reach the 
same growth (as) by today. Such efTects in linear growth 
and the CMB anisotropy power spectra allow that EDE 
must contribute less than a few percent to the early energy 
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density l|Doran. Robbers fc WetterichI [2007h . The acoustic 
sound horizon shifts by approximately half of this. However 
this article shows that observations of non-linear structure 
should not be expected to provide a substantial leap in our 
ability to detect EDE. 

This study has demonstrated t hat t he h alo mass func- 
tion fo rmulas of both lJenkins et al] (|200ll ) and lWarren et ahl 
l|2006h are valid for predicting the EDE to ACDM re- 
sults ratio to within ~ 10% or better in the mass range 
~ 10^^ - 10^^ Mrs,. Future studies of EDE should use these 
formulas or N-body simulations when considering halo abun- 
dances, rather than spherical collapse motivated mass func- 
tions as have been used in the past. 

We have als o demonstrated that the Halofit 
l|Smith et al.l l2003h formula for the non-linear power 
spectrum, when using the EDE linear matter power spec- 
trum and growth factor, is accurate for EDE cosmologies 
for scales larger than around k ~ 1/i/Mpc. On smaller 
scales this formula under predicts the power. Further work 
would be needed to accurately re-calibrate this formula for 
EDE on these s c ales, perhaps along the lines of the way 
iMcDonald et al.) (|2006l ) corrected Halofit for constant w 
dark energy cosmologies. Until such a study is performed 
we propose, as discussed in Section 14.41 that the results 
of this study suggest a simple prescription for predicting 
the non-linear power spectrum for EDE models to within 
an accuracy of a percent or so relative to the ACDM 
prediction. 
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